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Abstract. Let L = —A — W be a Schrodinger operator with a potential 
W S L^tt 1 (!>"), n > 2. We prove that there is no positive eigenvalue. 
The main tool is an L p — L p Carleman type estimate, which implies that 
eigenfunctions to positive eigenvalues must be compactly supported. The 
Carleman estimate builds on delicate dispersive estimates established in [7]. 
We also consider extensions of the result to variable coefficient operators 
with long range and short range potentials and gradient potentials. 



1. Introduction 



Let n > 2. Suppose W is a potential in M. n which decays at infinity. Then 
the Schrodinger operator 



has continuous spectrum [0, oo). In addition its spectrum may contain eigen- 
values which could be positive, negative of zero. Positive eigenvalues in the 
continuous spectrum are undesirable. They are very unstable since they are 
destroyed even by weak interactions between the continuous spectrum and 
the eigenvalue (see [9 ). Physically they correspond to trapped states in the 
continuous spectrum, and they are difficult to handle analytically. Moreover, 
excluding eigenvalues in the continuous spectrum is often a first step toward 
scattering. There is an extensive theory dealing with the absence of positive 
eigenvalues. 

It is well known that under weak assumptions like 



there are no positive eigenvalues. The argument uses Carleman estimates in 
three steps as follows. Suppose that 

—Aw - Wu = u 

with u e L 2 , where the eigenvalue is normalized to 1 by scaling. Then one 
proves that: 

(1) The eigenfunction u decays faster than polynomially at infinity. 

(2) If u vanishes faster than polynomially at infinity that u has compact 
support. 

(3) If u has compact support then it must vanish. 



W 



(1) 








These arguments work for many Schrodinger operators. However they do 
not cover Schrodinger operators for several particles (which are studied in j2] 
and PQ), neither do the standard arguments apply to the absence of bound 
states (i.e. L 2 solutions) in nonlinear optics modeled by problems of the type 

—Am = uju + a(x)\u\ a u 

with a bounded function a, because it is not clear how the assumption u G 
L 2 (R n ) is related to pointwise decay. 

On the other hand the assumption (Q) on pointwise decay is sharp: There 
is the famous Wigner-Von Neumann example of a positive eigenvalue and a 
potential decaying like l/\x\ but not better, see [121 El- 
Motivated by the above questions and by other potential applications one 
seeks to replace the pointwise bound (JTJ) by an LP bound. In terms of scaling 
any such bound must necessarily be weaker than ((TJ) due to counterexamples by 
Jerison and Ionescu (j^j) with potentials concentrated close to n— 1 dimensional 
planes. Jerison and Ionescu [3] have recently obtained absence of embedded 
eigenvalues for W G L n l 2 . In this paper we obtain the same result for a larger 
class of potentials which includes 

(2) WeL~ . 

We note that a higher index is better since it allows for potentials with less 
decay at infinity. Another way to look at this is that such a condition is mostly 
relevant for the low frequency part of W. The counterexample of Jerison and 
Ionescu (j3]) shows that this is the highest possible exponent. 

Our method is robust enough so that it also allows us to add a long range 
potential, and also to replace the Laplacian with a (mildly) asymptotically flat 
second order elliptic operator. The latter generalization is more technical and 
less self-contained, so it is discussed only in the last section. 

Thus we consider potentials which are the sum of weakly decaying long range 
potentials V and short range potentials W. We even include the eigenvalue 
A > into the long range potential and study the problem 

(3) (-A-V)u = Wu. 

To describe the long range potential we define the space C 2 X ^ by 

Definition 1. G\ L is the space of C 2 oc functions for which the following norm 
is finite: 

||/|U := max{sup \f(x)\, sup(x)|D/|, sup(x) 2 \D 2 f\} 

{x> X 

Then we introduce the condition 
Assumption A 1 (The long range potential). V belongs to Cf x v and satisfies 



(4) 



The bound from below on V corresponds to the condition A > while the 
last bound in says that for large \x\ the function \x\ 2 is strictly convex 
along the null Hamilton flow for —A — V, and thus guarantees nontrapping 
outside a compact set. 

To describe the short range potential we define the space 

1_ 2(n+l) 

Definition 2. X is the space ofW lo ™ +1 ' n+i functions for which the following 
norm is finite: 

\\W\\x = SUP ||W|X|| _ 1 2(n + l) /\\U\\ 1 2(n + l) U > 3 

U£C°° W n + ! ' n + 3 W n + l ' n — 1 

\\W\\ X = supJ\Wu\\ w _i +et e/\\u\\ w i_ ei6 n = 2,e>0 

For a domain D C lR n we denote 

X(D) = {1 D W; W G X} 

Then we introduce 

Assumption A 2 (The short range potential). W belongs to Xi oc and can be 
decomposed as W = W\ + W2 where 

(5) limsup ||W / i||x({xi2^<|x|<2^+i}) < 5 

j— t-OO 

(6) limsup |x||W2(V)| < 8. 

\x\^oo 

The W2 component corresponds to the L? Carleman estimates. The class of 
allowed W\ potentials includes Lt and or even better 1 £~(L?) where 

n + l 

the norm is taken with respect to a partition of lR n into unit cubes. 
Our main result is 

Theorem 3. Assume that V and W satisfy Assumptions Al and A2, let 
7~i > r and assume that 5 is sufficiently small. Let u £ H^JW 1 ) satisfy Q 
and (1 + |x| 2 ) T1 ~5u e L? . Then u = 0. 

By comparison, the result of Jerison and Ionescu applies to the case 
V = 1 and If 6 L2, n > 3. We note that the exponent p = n/2 is critical 
for weak unique continuation; for smaller exponents there are examples of 
compactly supported eigenfunctions, see 0. 

The conditions © and (jHI) have a different scaling behavior. Nevertheless 
both are sharp, which can be seen by the Wigner-Von Neumann example and 
the non radial counter example of Jerison and Ionescu. 

The proof uses Carleman estimates, following the same three steps indicated 
above. A combined L 2 - LP Carleman inequality replaces the previous L? Car- 
leman inequalities. Proving such inequalities is a highly nontrivial task and 

H§L 1+ if n = 2. 
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relies on the bounds established in [7j. Conjugation of the operator —A — V 
with the weight of the Carleman inequality leads to a non-selfadjoint partial 
differential equation. A pseudo-convexity type condition is satisfied, but it 
degenerates for large x. This is related to the fact that the anti-selfadjoint 
part of the conjugated operator decays for large x in relevant coordinates. 

Compared to earlier work and to the steps outlined above, we also consider 
a different family of weights in the Carleman estimates. Precisely, we begin 

with weights of the form h(x) = e T V / I^T for part El of the argument, which we 
then flatten at infinity for part ^ This yields a more robust argument, and 
also better results in the variable coefficient case. 

The paper is organized as follows. In the next section we state all the L p 
Carleman estimates and show how they lead to the result on the absence of 
the embedded eigenvalues. 

There are two main ingredients in the proof of the LP Carleman estimates. 
The first is the L 2 Carleman estimates, which are proved in Section El The 
second is a dispersive estimate for second order operators which is obtained in 
Section fusing an earlier result of the authors, namely Theorem 3 of [7]. This 
is of independent interest so we state it in more generality than needed here. 

The L p estimates are proved Section The L 2 bounds obtained earlier are 
used to localize the L p bounds to small spatial scales. Then we can rescale to 
a setting where the general dispersive estimates of Theorem |7] apply. 

Finally, in the last two sections we discuss the extension of the results to 
second order elliptic operators with variable but asymptotically flat coefficients 
as well as unbounded gradient potentials. This goes along the same lines. 



2. Carleman estimates and embedded eigenvalues 

As explained above the proof depends on Carleman inequalities. In this 
section we explain the Carleman inequalities and their application whereas 
most of the proofs are postponed to the remaining sections. 

Let 1 < p < oo and s6l We define the Sobolev space W s,p (R n ) by the 
norm \\f\\w* = 11(1 + \D\ 2 ) s/2 fhp and W S > P (U) for open subsets U of W 1 
through its norm which is the infimum of the norm of extensions. 

Given a measurable function / and the Sobolev space W s,q we define the 
norm 

i/p 



E 



p 

W s >i({23- 1 <\x\<2i+ 1 }) 



with the obvious modification for p = oo. 
Our Carleman estimates have the form 

i| e fc(bx([<B[)) v || i + iigMHM))™;!!^ < 

{ ) inf He^-IVVillw + lle^^/all x w 

/ 1+ / 2 =(--A-V> ™+a 



where p is given by 



min(|x|)) , ^(ln(NI)m(ln(|x|)) N^ 
\ |x| 2 |x| 4 

with /i" denoting the positive part of h". As a general rule, the function h is 
chosen to be 

(a) increasing, hi > r , with h'(0) large. 

(b) slowly varying on the unit scale, \h^'\ < h! for j = 2, 3, 4. 

(c) strictly convex for as long as /i'(ln(|x|)) > \x\. 

More precise choices are made later on for convenience, but the estimates 
are in effect true for all functions h satisfying the above conditions. 

The two terms in p have different origins. The second one simply measures 
the effect of the convexity of the function h. The first one, on the other hand, 
is due to the presence of the long range potential, which provides some extra 
pseudoconvexity for large \x\. 

A simplifying assumption consistent with the choices of weights in this paper 
is to strengthen (c) to 

(c)' h"(lxi(\x\)) « h'(ln(\x\)) for as long as h'(ln(\x\)) > \x\. 

This allows us to simplify the expression of p to 



(g) = ( h'(H\x\)) / + h'(ln(\x\)f 

\ \x\ 2 \ \x\ 2 

Our Carleman estimates use weights which grow exponentially, but also 
allow for the possibility of leveling off the weight for large enough \x\. 

Proposition 4. Suppose that V satisfies Assumption Al. There is a universal 
constant Eq such that with 

(10) h' e (t) = n + (rel - n 



t 2 + ee* 



(|7j) holds with h = h e for all \e\ < Eq, v supported in \x\ > 1 and satisfying 



\x 



T1 2v E L 2 , uniformly with respect to r large enough. 



The coefficient | in the exponent is chosen somewhat arbitrarily. However, 
it must be smaller than 1 in order for stage (c) above to be reached. This is 
necessary if we are to be able to taper off the weight at infinity. We continue 
with a short discussion of the weight h e . 

For small t it is uniformly convex in the sense that h" ~ h' e . The first 
interesting threshold for it is to defined by 



T 



This implies that h' e (to) ~ e to . In the range [0, to] the last factor in fllOl) 
is largely irrelevant, and h' e behaves like an exponential. In this region, the 
pseudoconvexity in the Carleman estimates is produced by the convexity of h. 



After to h t is still convex, roughly up to t\ defined by 

e* 1 « e _1 r 2 

The region t\ + 0(1) contains both the inflexion point t\ and the maximum 
point for h' e . In between to and ti the pseudo-convexity comes from the poten- 
tial term, while the contribution from the convexity of h e is still positive but 
smaller. 

Beyond t\ + 0(1) the function h' e (t) — t\ decays in an exponential fashion. 
The last interesting threshold is £2 where h' e approaches 1, given by 

Between t\ and ti there is still convexity coming from the potential V, which 
suffices in order to control the lack of convexity of h e . Finally, after ti the 
pseudoconvexity in the classical sense is lost, but there remains an Airy type 
gain to push the estimates through. 

Proof of Theorem^ Here we show that Proposition 0] implies Theorem |H| 

STEP 1: We prove that u decays at infinity faster than e~ T v^. We choose 
R large enough so that (see Assumption A2) 

(11) SUp ||Wi|| X({:r |2 J <| 2: |<2 J +i}) < 25 

2J+ 1 >_R 

(12) sup |x||W 2 (x)| < 25. 

\x\>R 

Choose cf) G C°° be identically 1 for \x\> 2R and for \x\ < R. We set v = <ftu. 
Then 

-Av - Vv = Wv - (A(f))u - 2V0 • Vu 

For t\ as in Theorem El we have |x| T1_ ^f G L 2 , therefore we can apply Propo- 
sition 0] with e > to v to obtain 

11^(1*1*1)^1 +1^^1-1)^11^ < \\e h ^Wp-\\u\ + \Vu\)\\ lHB2rXBr) 

+ \\ e h*(}n\ X \) w u + \\ e h4^\x\)-l W U 

By (fTT|) . (|T2*|) if 5 is small enough then we can absorb the last two right hand 
side terms on the left to obtain 

\\e Mln]x]) v\\,r Xem + II^^Hb S H^'V'tM + \Vu\)\\ lHB2R \ Br) 

Then letting e — > in the definition of h yields 
(13) 

HeVMuH ^ + \\e^pv\\ L 2 < \\e T V^p-\\ u \ + \Vu\)\\ lHB2R \ Br) . 

which shows that v and therefore w is rapidly decaying at infinity. 
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STEP 2: We prove that u vanishes outside a compact set. This is done 
using (|13|) (which can also be derived directly from Proposition 0] as above). 
From ()13j) we obtain 

JrV^||eVWi;|| _ + ||e^V|b < ||(| w | + | Viil)^^^. 

Letting r — > oo shows that t> = outside _B 2 i?- Then the same holds for u. 

STEP 3: We prove that u is identically 0. Assume by contradiction that 
this is not the case, and choose r minimal so that u is supported in 5(0, r). 
Our problem is scale invariant, so without any restriction in generality we 
can assume that r > 1. Take Xq G supp u with \x§\ = r. The problem 
is also invariant with respect to translations so we can assume instead that 
supp u G B(x ,r) and 2x G supp u. 

To reach a contradiction we prove that there is a > so that u is supported 
in B(0,2r — a). This follows as in STEP 2 provided we know that for every 
5 > we can find p > such that 

W w ^ v \\ w -^m - 5 iHiv^^' supp v c B{2x ^p) 

Then a is chosen so that 

{2r - a < \x\ < 2r} n 5(a;o, r) C B(2x , p) 

Due to our choice of W this is a somewhat technical matter which is left for 
Proposition El in the appendix. This step can be approached alternatively by 
the unique continuation results of [Zj. □ 

3. The L 2 Carleman estimates 
In this section we obtain the L 2 Carleman inequalities. 

Proposition 5. Suppose that V satisfies Assumption Al. Let h be as in ffTTTf) 

and p as in (jHJ). Then for all u satisfying |a;| T1_ 2u G L 2 we have 
(14) 



|e h(ln|a VlU 2 + 





\x 




h'(ln | 


x\) + 


x\ 



<|| e MlnN) p -l( A + y )u || i2 . 

L 2 



uniformly with respect to r sufficiently large and < e < 6q. 
Proof. We use a conformal change of coordinates 

t = In \x\, y = x/\x\ G § rt_1 

Denote 

Au = g 

and set 

v(t, y) = e^Me'y), f(t, y) = e^^g^y) 
A routine computation shows that 

N (n + 2)/2 (A + V)lx \(n + 2)/2 = <?_ + ^ _ {{n _ 2)/2) 2 
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therefore v solves the equation 

(15) Lv = f, L = d 2 + As»-i - ((n - 2)/2) 2 + e 2t V 

We also note that part of Assumption Al in the new coordinates we get 

. , V t 1 
— hm ml — = r < - 

t^oc 4V 2 

By (J1J) we slightly readjust r and choose t so that 

(16) < r < -, t > t 
For any exponential weight h we have 

(17) [ e 2h ^\u\ 2 dx = [ [ e 2fc W- Hrt |tt(t 2 /)| 2 dtdj/=||e h We t t;||i a(Rx5 „_ 1) , 



e 



2 ^I-D|^| 2 dx= / / e 2 ^b(iy)| 2 ^dy=|| e ^e-V||l 2(RxS .- 1) . 



Hence, in the new coordinates the bound ()14|) becomes 

(19) \\e h U pi v\\v + ||e fe W ef+ ^ /(t) V^|| L2 < \\e h ® p^fU*, 

where Vf is the gradient of v with respect to y and t and, by 0, 

pi (t)=e t p = h'{t) (e 2t + ti(t) 2 )* 

To prove the above bound one would like to follow a standard strategy. This 
means conjugating the operator with respect to the exponential weight, and 
producing a commutator estimate for the self-adjoint and the skew-adjoint part 
of the conjugated operator. There are two small problems with this approach, 
both of which occur in the region where h'(t) is small. 

First we want to incorporate the weight p^ 1 on the right, which would 
require an additional conjugation. Where h' is small this cannot be treated 
as a small perturbation, so we really have to include p~ l in the exponential 
weight. 

This leads to a second difficulty After including p _1 in the exponential 
weight the commutator between the self-adjoint and the skew-adjoint part 
of the conjugated operator is no longer fully positive definite and we need a 
slightly modified argument. 

To handle both issues we prove a slightly more general result and then we 
obtain (J19)) as a special case of it. Precisely, we consider an exponential weight 
<p as follows: 

(i) (j)' > Tx — |, and (f)'(0) is large. 

(ii) 1 + 0' is slowly varying on the unit scale, i.e. 

|<^)(t)|<l + 0'(t) j = 2, 3 

(iii) 4>' can only have a limited exponential growth rate, 0" < 4(1 + (j)'). 
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Together with (i) this yields the existence of a unique t so that (f>'(t ) = e*°. 
Our last assumption asks for uniform convexity up to to : 

(iv) 4>"{t) ks 4>'(t) for < t < t + C for some large parameter C. 
We summarize the bound for the weight e^: 

Lemma 6. Consider a weight function satisfying the conditions (i)-(iv) 
above. Then for all v which are supported in t > and with e^ +t v G L 2 we 
have 

(20) \\e^ l \e 2t + <p\tf)^v\\ L 2 + ||e^V^|| L 2 < ||e^(l + 0')~*^IU»- 



Proof. First we conjugate with respect to the exponential weight. If we set 
w = e^v then w solves the equation 

We decompose Lh into a selfadjoint and a skewadjoint part, 

1$ = % + A - (^) 2 + e 2 '\/ + 0' 2 , LJ = -0'9 t - 9,0' 
The bound to prove is 

(21) IK^ + ^O^'wH^ + UVHIl" < IKi + ^T^^IU- 

The proof of this inequality is based on several integrations by parts. In a 
standard manner one verifies that the integrations by parts below are valid if 

e <l>+t v e L 2 

We multiply P^w by — \w t and integrate by parts to obtain 

/ <f>'\w t \ 2 dydt+ f{-S'" + \<t)'(t>")\w\ 2 dtdy+ [ ^-(2V + V t )w 2 dy dt 

(22) J If 

= - wtL^wdydt 

This computation is essentially like taking the commutator of and L\. On 
the left we have mostly positive contributions, with the following qualifications: 
-the first term can be negative where 0' < 

-the <f>'<f>" term can also be negative, but only for t > t + C where it is 
controlled by the V term. 

-the (j)'" term is controlled either by the V term or by the (fJ(f>" term. 

To correct the first term in the region where 0' is negative we consider a 
cutoff function x which equals 5 in {<fi' > 2} and which equals 1 in {0' < 1}. 
Here 5 is a small universal parameter which we shall choose below. Since 0' + 1 
is slowly varying we can assume that x has uniformly bounded derivatives. 
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Multiplying P^w by x 2 (t) w an d integrating gives 
(23) 



X^lli 2 + llxV^||i 2 + (^) 2 ||^||i 2 - / x 2 {e 2t V + <P' 2 )\w\*dydt 



-{d 2 x 2 )w 2 dy dt + J wL^w dy dt. 

We multiply this by \i and add to the previous relation. This yields 
(24) 

A*||xVtt>||£ 3 + J(x 2 V + <f>')\wt\ 2 dydt+ J Q-XV + ^J e 2t Vw 2 dydt 

+ l^'4>"-x 2 ^' 2 )H 2 dydt 

[ - 0'74 - (^) W + fl£xV] M 2 dy * 

+ y \x 2 ^w + ^w t )L (j> wdydt. 

To ensure that the left hand side is positive definite we recall that for large t 

Vt 1 „ 

-W- T0<Tl - 2 +<P 
Hence if we choose /i positive so that 

1 1 

2-n</i<--r 

then the first three terms are positive definite. 

For the fourth term we consider two possibilities. If t < to + C then x — $ 
while <fi" ~ cf)' so it yields a positive contribution. We choose the universal 
constant 5 so that 

if t < t + C. For larger t this fourth integrand may be negative but then it is 
controlled by the third. The first term on the right hand side is controlled by 
the left hand side and we obtain 

||VH|!a + IK 1 + 4>')*™tfv + \\(<f>'(t) 2 + e 2t )^w\\ 2 L2 < J(fxw + \w t )L^wdydt 

The proof is completed by an application of the Cauchy-Schwarz to the right 
hand side. 

□ 

Proof of Proposition continued. 

We obtain (JT§j) from Lemma El For this we need to associate to each weight 
h a function satisfying (i)-(iv) with the property that 

1 + (j)' « ti, (1 + 0')~*e* « e h (h' 2 + e 2t )-^ 
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The natural choice for is 

0(i) = h(t) -t + l ln(l + h'it)) - \ ln(l + e-*h'{t)) 

Then 

1 h" {h' - /i")e"* 



0' = /i' - 



2 4(1 + h') 4(1 + e-*/i') 

We verify the properties of 0. It is easy to see that 1 + 0' is slowly varying. 
This implies that the last two terms in 0' are bounded and have bounded 
derivatives. Hence the properties (ii)-(iv) follow from the similar properties of 
//'. 

It remains to check the bound 0' > T\ — |. This is clear when b! 1 which 
corresponds to ee^ <r 3 . For larger t we have 

h'(t) = r 1 + K 3 e-Hl + 0(r- 1 )) 

and 

/ i "(t) = -lrVl(l + 0(r- 1 )) 

Then 

0'(t)>r 1 -i + irVl(l + O(r- 1 )) 

so the desired bound is again verified. We note that what happens when 
h! is small is not so important anyway; in this region we can simply choose 
0(t) = h(t) □ 

4. A GENERAL DISPERSIVE ESTIMATE FOR SECOND ORDER OPERATORS 

In this section we study the second order operator 2 

L M = d i a ij (x)d j + fi 2 c(x) - i^{b 3 {x)d j + djV{x)\ 

in the unit ball B C M ra , n > 2 with real coefficients a 1 - 7 and complex coefficients 
ti> and c. Here /i is sufficiently large and plays the role of a semiclassical 
parameter. Concerning the type and regularity of the coefficients we assume 
that 



(REG) 



the matrix (a y (x)) is real, symmetric and positive definite 
the functions a iJ , b % and c are of class C 2 



We define the symbol 

= C""t.r)i ; + c(x) + 2bjZj 
The real part of / is a second degree polynomial in £ with characteristic set 
charM(x,0 = {£ G R n ; 3W(a;,0 = 0} 



We use the summation convention here and in the sequel. 
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The geometric assumption on the operator L is 

, , J for each x the characteristic set char x Ml(x,^) 

^ ' 1 is an ellipsoid of size ~ 1. 

Our third hypothesis is concerned with the size of the Poisson bracket of 
the real and imaginary part of L. We are interested in a principal normality 
type condition of the form 

(25) |{3M(z,0,9tf(z,0}l;$* + l^foOl + 

where the relevant range for 5 is pT 1 < 5 <C 1. This would suffice for our 
purposes if in addition we knew that all the coefficients of / are of class C 3 . In 
general for technical reasons we need to replace the inequality with a decom- 
position 

(26) $Sl}(x, = Sq (x, + q{(x, 03W(*, + q\{x, £)$Sl(x, + q 2 (x, £) 
Thus our last assumption has the form 

the Poisson bracket ^sl} admits a representation (|2fjj) where 
(PN) { \d«dP qi (x,0\<c al 3 \a\<i 

ko|<i, |gii + |gil<i, |g 2 |<|/| 

For L in the class of operators described above we are interested in con- 
structing a parametrix T which has good LP — > L p and L 2 — > L p mapping 
properties, while the errors are always measured in L 2 . A dual form of this 
also allows us to estimate the L p norm of a function u in terms of the L 2 norms 
of u and Lu. 

In the context of the Carleman estimates such parametrices allow us to 
superimpose local L p — > LP bounds on top of the global L? — > L? estimates 
in order to obtain a global LP — > LP bound. 

Such estimates are dispersive in nature and are strongly related to the 
spreading of singularities in the parametrix T. This in turn is determined 
by the nonvanishing curvatures of the characteristic set char x $ll(x, £). 

If L has constant coefficients and real symbol then the theorem below is 
nothing but a reformulation of the restriction theorem. If L has real symbol 
but variable coefficients then we are close to the spectral projection estimates 
of C. Sogge In the case when L has constant coefficients but complex 

symbol some bounds of this type were obtained in 0j. 

In the more general case considered here we rely on bounds and parametrix 
constructions in the author's earlier paper jjj. These apply to principally 
normal operators. The operator L M is principally normal on the unit spatial 
scale only if 5 ~ pT 1 . Otherwise, we use a better spatial localization to the 
(5/x)~2 scale. On one hand is principally normal on this scale, while on 
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the other hand this localization is compatible with the L 2 estimates and this 
allows us to easily put the pieces back together. 

All Sobolev norms in the theorem below are flattened at frequency \x instead 
of frequency 1 as usual. Hence we introduce the notation 

W* = {ue S'; (fi 2 + D 2 )iu G LP} 

with the corresponding norm. 

We note that the operator L is elliptic at frequencies larger than fi so all 
the estimates are trivial in that case. All the interesting action takes place at 
frequency < /i, where we can identify all Sobolev norms with L p norms. 

Theorem 7. Suppose that the operator satisfies the conditions (REG), 
( GEOM) and (PN) for some 5 > ^T 1 . Let <p G C(i?2(0)) have compact support. 
Then 

A ) There exists an operator T such that 

\\Tf\\ i ^ + (ty) 1 V 1/2 ||T/||*i 

(27) 

< . mf (^r V V 1/2 ||/i|U 2 + II/2II 1 »fe±» 

f=fl+J2 w„ n+a 

I- 1 

and 



(28) mf (^-^V^IIAII^ + IIMI 1 afew 

/=/l+/2 Wfj, ™ TT ' n+ ' d 

B) For all functions u in B 2 {Q) we have 



\\M\ 1 w) <(^) x V /2 lkl 



(29) 



n + l ' 71 — 1 



+ inf f^-'V^II/illr' + ll/all i ^ 

Lu=f 1 +f 2 w n+l w+3 



Suppose that in addition the problem is pseudoconvex in the sense that 
(30) Qb(ar, ~ 5 > /i -1 x G £ 2 (0), /i > 1 

T/ien /or all functions u with compact support in -82(0) we have 

\\U\\ 1 2(n+l) + {b~H) 1,A H l/2 \\u\\ L 2 

(31) 

< mf (M^V^II/iIIl^ + II^II 

Lu=f 1 +f 2 +f3 w "+!' «+3 

The difficult part of this theorem is the existence of the rough parametrix in 
Part A. This existence will be derived from Theorem 3 in [7j. The arguments 
repeat partially those of Section 3, 7 and 8 of jjj. 

Proof. Part A. (i) Localization. We first reduce the problem to the case 
when 8 = This is done by localization to a small spatial scale and then by 
rescaling. The appropriate spatial scale is r = (//<*)-*. We cover the support 
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of (f) with balls Bj of radius r and choose a subordinate partition of unity of 
the form 

Suppose that within Bj there exists a parametrix Tj satisfying the desired 
estimates. Then we set 

N 
3=1 

The bound (|27j) for T follows directly by square summing the similar bounds 
for Tj. For (J2Hj) we compute 

N N 
3=1 3=1 

For the first term we use (}2*%j) for Tj while for the second we estimate the 
commutators using (|2*7)l for Tj. 

In order to obtain the localized parametrices Tj we rescale Bj to the unit 
scale. Then the problem reduces to the original one but with 5 = fi^ 1 . 

(ii) The elliptic high frequency parametrix. 

For each x the zero set of 3?/ is an ellipse contained in a ball of radius -Br m (0) 
with R ~ 1. Let ip G C°°(IR") be a nonnegative radial radially decreasing 
function supported in B 2 (0) and identically 1 in -Bi(O). Let be as in the 
statement of the theorem. We fix a nonnegative function <fio G C°°(i?2(0)), 
identically 1 on the support of 0. We define Thigh by its Weyl symbol 

Then the following I? bounds are immediate: 



\\Thighf\\m 



< 



H7 t 



\\(l-LT high )(l-i>(D/(2 f MR)))<f>f\\ L , < H/ll^ 

This estimates are the elliptic versions of the parametrix bounds. By Sobolev 
embeddings they imply bounds of the type of Theorem 4. 

(iii) The low frequency parametrix. We first mollify the coefficients of 
on a scale fi" 1 ' 2 and note that this does not affect the hypothesis of the 
Theorem. We also modify its symbol for large £ and extend it to M. 2n so that 
it is of size /i 2 and so that it satisfies 

\O x 0^f,{X,l;)\ £ | ^i+H/2-|/3| jf j a j > 3 
By Theorem 3 of [7] there exists a parametrix Ti ow for Z M satisfying 

fi^WT^fW 2(n + l) + H 1/2 1 1 T, nw f I \ L 2 

(32) L ^ TT ~ 

< inf u- 1 / 2 ||/ 1 || L2+/ ,-^i||/ 2 || ^ 

f=fl+j2 L n=l 
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and the error estimate 

f ,- 1 / 2 \\(l-i;(x,D)T low )^(D/(2^R))<j ) f\\ L , 

m < inf ^WhWv +vtMM\ t 

f=h+h L s=i 

(iv) The complete parametrix In the final step we combine the low and 
high frequency parametrices. We set 

T = T high (l - + MiD/^Tiow^D^fiR)^ 

The estimate (|27j) follows easily from the similar bounds for T high and T iow . It 
remains to consider the error estimate. We have 

{I-LT)4>f = {I-LT high ){l-iP{D/2nR)(t>f 

+ ^(D/AfiR^I - L^T^iD^R^f 

+ [L»,M(D/4pR)]T low i/;(D/2tiR)<l>f 

+ (L- L™)M(D/tyR)T low il>(.D/2tiR)<j>f 

For the first two terms we use the error estimates for T high , respectively Ti ow . 
In the third term the commutator has size /i in L 2 so we can use the L 2 bound 
for T high . The operator 

(l - l;)m(d/^r) 

also has size \i in L 2 since the original coefficients differ from the mollified ones 
by /x -1 . This complete the proof of the inequality (|2Bj). 

1_ 2(n + l) 

Part B. We prove flU by duality as in Section 3 of [7j. Let <? G n+1 ' n+3 . 

We decompose <pg as 

0# = h + L*T(f)g 

where T is the operator of Theorem constructed for the formal adjoint op- 
erator L*. By part A of the theorem we have 

(^V^IWb + (^) 1/4 /i 1/2 ||T^|| i2 + ||T^|| 

< HSU _ 1 2(n+l). 

Therefore we can write 

l(<M,#)l =\{u,<f>g)\ 

<\(u,h)\ + \(u,L*T(Pg)\ 
=\{u,h)\ + \{Lu,T(f>g)\ 

<f(^) 1/ V /a ||«llL»+ inf (^-'V^ll/illxsi+ll/all . 

\ /J, 
X ||p|| 1 2(n + l) ■ 



1 2(n + l) 



n+3 



This implies the estimate 
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Part C. We begin with an L 2 estimate. The principal symbol of 

L, = + \Df)^ 

is 

Ux, = ( - ^(xMj + V 2 W(x) + 2/V&) (p 2 + |£|T 1/2 - 
A short calculation shows that 

5fi - {m^x, o, $*Ux, £)} < \Ux, o | 

and hence, by Corollary 11.14 of [TT], we obtain the bound 

<MIHU 2 ~ II^HU 2 + IMU 2 

If 5/i 3> 1 then the norm of u on the right hand side can be hidden on the left 
hand side. Applying this to w = (/i 2 + \D\ 2 ) l l 2 v we obtain 

(34) ^\\v\\ 2 < || V,||| 2 

For u as in the theorem we write 

u = v + TL^u 

The bounds for the second term come from part A. On the other hand, 

L^v = (1 - L^L^u 
for which we can use the error estimate (l28l) to obtain 



{5^y x '^- x ' 2 \\Lv\\ L ,< mf f^-'V^ii/iii^ + iiMi 

Then we successively apply (p£ij) and (|2*U|) to t>, concluding the proof. □ 

5. The L p Carleman inequality 

In this section we prove Proposition 0] We first conjugate with respect to 
the exponential weight. If we set w = e h ^ ln ^ x ^v then we can rewrite (J7|) in the 
form 

|HI i 2 ( n+i) +\\pw\\ L 2 < inf ||p _1 /i||l 2 + 11/2 || 1 3(n+l) 

£2^^' n-i L h w=fi+f 2 l 2 W~^' ™+a 

where 

L h = A + y W + //(In |x|) 2 |x|- 2 - /i'(ln |x|) [V^r + ^Vl 

We want to apply Theorem on dyadic annuli 

Aj = {x\2 j - 1 < \x\ < 2 J+1 } 
The rescaling y = 2~- 7 x transforms this set to Aq and the operator Lh to 
L{ = A + 2 2 >V + ti(H2>\y\)) 2 \y\- 2 - h!{\n{V\y\)) [V^ + JLy] 
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We verify that we can apply Theorem [7| to L? h . Since hi varies slowly on the 
unit scale we can take the corresponding value for p to be 

The coefficients b and c are given by 

c = tf(2*v + h'(H2'\y\)?/\vn t, = J^MlJt 

to \y\ 

and are clearly of class C 2 and size 0(1). We have 

Their Poisson bracket has the form 

,2 , . , tl h'(t), ltl2 , , , V 1 fc"(t) 



\y\ 2 Vj MVj 

Then we can apply Theorem [7| with 5 comparable to the size of the third term. 
For our choice of h we have \h"\ < h' and also 

h"(t) < < e* 

Hence we can choose 

5, = pf (2 2 ^'(jln2) + / i / (jln2) 2 / i / ;(jln2)) 
Let e Cq°(M) be a nonnegative function supported in [—1, 1] with 



J=-0O 



and let = 0(ln |x| — j). After rescaling, part A of Theorem [7| yields a 

parametrix Tj for in Aj with the property that 





\x\ 




h'(\n\ 


x\) + \ 


x\ 



l T ^t^%^ + l|pT ^ lli2 + ]lp h'an\ x \) + \x\ v{TML2 



+ \\p 1 (L h T j -l)4>jg\\ L 2 < inf \\p VlIU 2 (a,) + \\92W i a(n+i) 

9=91+92 W n + l' n+3 (Aj) 

We define a parametrix for Lh by 

oo 
3=0 

Summing up the bounds on Tj we obtain a bound for T, 

\\ T 9\\^i wi + \\pTg\\ L 2 < inf \\p' l gi\\ L 2 + ||# 2 || i w> • 
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We also compute the error 

oo oo 

1 - L h T = ^ 'M 1 - L hTj)4>j - ^2[L h , (pj]Tj4>j 

j=0 j=0 

Since 

[L h , <P,} = 0{\x\- l )V + 0(h'(\n \x\)\x\- 2 ) 

and 

ai (m |a;|j + \x\ 
we can bound the error by 

\\p-\l - LT)g\\ L 2 < inf Hp" 1 ^!!^ + H&H i a(«+i) ■ 

9=91+92 l 2 W nTT' "+ a 

Now, after the construction of the parametrix the assertion of Proposition 
0] follows exactly as the corresponding part of Theorem [7| We repeat the 
argument. Split w into 

w = v + TLw 

Then the second term satisfies the desired bounds while for the first we know 
that 

\\p~ l Lv\\ L 2 = \\p~^{LT - || £a < inf H^VIIl^ + ||<b|| i a(n+p ■ 

Lw=gi+g2 l 2 W ™+a 

Lemma El allows us to also estimate 

IIHU 2 

On the other hand by Theorem [?l B rescaled and applied to v in Aj we get 



i 5, ||HU 2 (a.) + Hp 1 Lv\\ L 2( Ai) 

W n+1 ' n— 1 



and after summation in j 
thereby concluding the proof. 



*>ll i < IIHU 2 + ||p || £ a 



6. Equations with gradient potentials 

In this section we discuss the corresponding results which are obtained when 
short range gradient potentials are added. Thus we consider equations of the 
form 

(35) (-A - V)u = Wu + Z l Vu + VZ r u 

with V and W as before. The gradient potential Z = (Z l ,Z r ) is subject to 
the following conditions: 
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Assumption A 3 (The short range gradient potential). The gradient potential 
Z G l°°{L n ) satisfies 



In addition for some R 3> ||V||.z,oo the low frequency part S < rZ of Z satisfies 
the conditions in Assumption A2. 

The L n assumption is natural due to scaling. The low frequency condition 
is also natural, since on the characteristic set of —A — V the frequency has 
size 0(1), and at frequency one there is no difference between the potential 
and the gradient potential. Under these conditions we have 

Theorem 8. Assume that V , W and Z satisfy Assumptions A1,A2 respec- 
tively A3. Let 7~i > To and assume that 5 is sufficiently small. Let u G H^JW 1 ) 
satisfy © and (1 + \x\ 2 ) T1 ^u G L? . Then u = 0. 

By scaling we obtain the following result on the absence of embedded eigen- 
values: 

Corollary 9. Assume that V , W and Z satisfy Assumptions A1,A2 respec- 
tively A3 with 5 = 0. Then there are no embedded eigenvalues for the operator 



The problem of introducing gradient potentials has long been considered 
in the context of the unique continuation and the strong unique continuation 
problems for the same operators as here. There the key breakthrough came in 
Wolff's work 13] who proved that Z G L n suffices for the unique continuation 
property. He also obtained the same result for strong unique continuation 
but only in low dimension. Later his ideas were used by the authors in j3] to 
complete the picture for strong unique continuation in high dimension, working 
with gradient potentials Z G l x L n . This latter paper is more relevant to the 
present context as it provides Carleman estimates in largely the same format 
as here. 

Ideally, one would like to include matching gradient estimates to our L p 
Carleman inequalities. This would solve the problem but unfortunately cannot 
work. Wolff's contribution was to show that by osculating the weight one can 
considerably improve the bounds for the gradient term in the equation. Thus 
the choice of weights ultimately depends both on the gradient potentials and 
on the solution u. In our context this argument is needed only at spatial scales 
where the frequency of the conjugated operator is larger than one. Elsewhere 
the gradient does not contribute much to the problem. Thus we are led to 
consider perturbed weights 



(36) 



limsup ||Z||zn({ a .|2j<[ !r |<2j+i}) < 5 



A - W - Z l V - VZ r 



(37) 



i[) etT (x) = /i e (ln \x\) + k(x) 
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where k is not spherically symmetric but is small in an appropriate sense. The 
assumptions on k are summarized in what follows: 



(3* 



supp k C {\x\ < r 2 } 

\x\ a \V a k(x)\ < h'{ln\x\) a = 1,2, 3 



Part of the Carleman estimates below describes what happens in elliptic 
regions of the conjugated operator L^. To select (part of) this elliptic region 
we introduce a pseudodifferential operator x>r which selects the region 

E = {\x\>t\ \i\>R) 

Here both the truncation in x and in £ are done on the dyadic scale, while R 
is chosen sufficiently large so that E is away from the characteristic set of P^. 
Then the Carleman estimates are as follows: 

Theorem 10. Assume that the long range potential V satisfies Al. Let Z 
satisfy A3 with 

H^H^l™ + \\S<rZ\\i< x x+(x)- 1 l< x < 1 
Then for each < e < eo, t large enough, T\ > r and v which satisfies 
(1 + |:r| 2 )' ri ~2t! g L 2 there is a weight perturbation k satisfying (J37|) . (|3*%|) so 
that the following estimate holds with constants independent of < e < eo, 
r > r : 

(39) 

He^Ml -J- *c» +1 ) + \\x>Re^ A Mm + \\e^ x) pv\\ L , 
+ \\e^ {x) Z l S7v\\ x 2( n+i) + \\e^^ x) VZ r v\\ , M 

PW "+T' n+3 +X>R H~ 1 l 2 W n + 1' « + 3 +X>R H~1 

< inf lle^rWp-i/.ll^ + lle^rW^n x 

The key feature of the theorem is that the weight ip e)T (x) depends both on 
the potential Z and on the solution v itself. Once this result is established, 
it leads as before to the conclusion that solutions to (|33|) must be compactly 
supported. Then (a variation of) Wolff's weak unique continuation result ^3] 
takes over and implies that v must be identically 0. We also refer the reader 
to [3] , where the estimates are formulated in a way similar to this paper, and 
where both left and right gradient potentials are considered. 

The proof requires the following steps: 

(i) Conjugate the equation with respect to the exponential weight and set 
w = e^ e ' T v. This elliminates the exponential weight from the equation and 
replaces the operator —A — V by its conjugated operator L^ T . 

(ii) prove the L 2 estimate for v uniformly for all weights <p ttT with k satisfying 
(JHHj) . This is done exactly as in Section El The size of the perturbation k is so 
that its effect is negligible in this computation. 
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(iii) prove the estimate (J7j). again uniformly with respect to all choices for 
k. This repeats the arguments in Section with no change. 

(iv) Add in the H~ l and H 1 norms, thus proving (J3H|) for Z = 0. This is 
done in an elliptic fashion, by constructing an elliptic parametrix for L$ away 
from its characteristic set. For this the norms involving p are used only to 
estimate errors, while the LP norms are all used via L 2 norms and Sobolev 
embedding. The standard pseudodifferential calculus can be applied since the 
coefficients of are smooth on the dyadic scale in x. 

(v) observe that the L 2 estimates allow localization on the dyadic spatial 
scale. Thus we separate the estimate into two regions, {|x| < r 2 } and {\x\ > 
r 2 /2}. 

(iv) show that within the first region it is possible to choose the weight 
k so that the estimate with Z included holds. This is the part that uses 
Wolff's osculation lemma, and it is explained in detail in 0. Our case here is 
somewhat simpler than in [S] since in the region {\x\ < r 2 } we have uniform 
convexity of the weight, h" w h! . Also the LP bound here is stronger than in 

which only makes things better. 

(v) prove the estimate in each dyadic component of the second region \\x\ > 
r 2 /2} with Z included. This starts from the estimates without Z and uses only 
elliptic bounds. We outline the argument. Since we use dual norms on the left 
and on the right of (J39|) . it suffices to do it for the ZV term. The bounds for 
VZ will work out similarly but in dual spaces. 

We split Z into a low and a high frequency part, 

Z = Z <R + Z >R 

and the gradient also, 

V = V<R/2 + V >J R/ 2 

77.4-1 

Using the bound on Z <R / 2 we can directly estimate the contribution 
of Z <r V < r/2 which is located at low frequency. 

The contribution of Z >R V <R /2 lies at high frequency, so it suffices to bound 
it in H~ x . We can actually bound it in L 2 , 

\\Z>rV <R / 2 W\\ L 2 < \\Z\\ L n\\w\\ 1 2(n + l) 

n-1 

For ZV >R /2 we can use the H 1 bound to write 

\\ZV >R/2 W\\ 2n < \\Z\\ L n\\x >R /2W\\ H l 

and conclude by Sobolev embeddings. 

7. Asymptotically flat metrics 

In this section we describe how the results on the absence of embedded eigen- 
values extend to variable coefficient asymptotically flat metrics. We replace 
the Laplacian with a second order elliptic selfadjoint operator 

L = ~d ja jk d k + i{Vdj + djV) + c 
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where the coefficients a, b, c are real. We assume that P is flat at infinity in 
the sense that (see Definition [TJ : 

limsup I^UVa 1 - 7 ! < 5o, 

(40) |x|^oo 

limsup \b(x)\ + \x\\Vb(x)\ < 5\, limsup \c(x)\ + |x||Vc(x)| < 5 2 

\x\— >oo |x|— >oo 

We also slightly strenghten the assumption A3 to make it stable with respect 
to changes of variable: 

Assumption A 4 (The short range gradient potential). The gradient potential 
Z G l°°{L n ) satisfies 

(41) limsup \\Z\\ Ln{{x \ 2] <\ x \< 2] +i }) < $ 

In addition (D)~ N Z satisfies the conditions in Assumption A2 for some N 
sufficiently large. 

Then we have 

Theorem 11. Assume that W , V and Z satisfy Assumptions A1,A2 and A4 
with small enough 5, that T\ > r and that the coefficients of P satisfy (|4U|) 
with So and 5± sufficiently small. If u G i// oc (lR n ) solves 

(42) Lu + Vu = Wu + Z l Vu + VZ r u 

and (1 + \x\ 2 ) Tl ~^u G L 2 then u = 0. 

The assumption of Theorem ^2 are not scale invariant. For the following 
straightforward consequence we rescale the operator. 

Corollary 12. Assume that the coefficients of the operator P satisfy ()40|) with 
5q sufficiently small. Let W , Z be as in Assumptions A2, A4 with 5 = 0. Then 
there exists C > so that P + W has no eigenvalues A > C8\. 

The proof follows the same outline as in the constant coefficient case. We 
describe the steps in what follows, and discuss the necessary modifications. 
First one needs to augment (J4(Jj) to gain also the relation 

(43) limsup \a(x) — I n \ < 5 

\x\— >oo 

This is achieved using a change of coordinates somewhat similar to the one 
introduced in Due to (J4T?|) . within each spatial dyadic region this can be 
achieved with a linear change of coordinates. But from one dyadic region to 
the next these linear maps differ by 0(6). Hence gluing them together yields 
a nonlinear function \ which achieves ()43j) and has the regularity 



l^xO*)! <<5iM 1_H l«l>2 
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It is easy to verify that such a change of coordinates does not affect 8\ by more 
than a fixed factor. 

If x were linear then the Assumption Al on V would rest unchanged. As it 
is, we have to modify t\ by 0(8q), which is suitably small. 

Finally, the operator L is still L 2 selfadjoint in the new coordinates but with 
respect to the measure given by the Jacobian J of the change of coordinates. 
This implies that JL is selfadjoint with respect to the Lebesque measure. This 
requires replacing V and W by JV and JW, which has no significant effect 
on our assumptions. 

Once (J43j) is gained the Carleman estimates (J7J) in Proposition 0] remain 
valid with essentially no change. The only minor modification that is needed is 
concerned with what happens within a compact set, where we have no control 
over the geometry of the coefficients a y ' in the principal part. But this can 
be easily addressed by adding some additional convexity to the exponential 
weight within this compact set. Precisely, a weight of the form 

h{t) = re xt 

would suffice for bounded t provided A is large enough. 

The L? Carleman estimates are established using integration by parts, and 
do not require any bounds on the second derivatives of a iJ . 

The LP Carleman estimates are derived from the L 2 ones exactly as in 
Section For comparison purposes, we recall that the LP estimates proved in 
[7j and [S] only require bounds on the first derivatives of a l K This is because 
the spatial localization which is allowed by the Carleman estimates is on a 
scale on which one is allowed to freeze the a*- 7 modulo negligible errors. The 
same applies here for \x\ < r 2 (which corresponds to e* < r 2 ). However, 
beyond this threshold the rescaled skewadjoint part becomes very small and the 
problem is close to the spectral projection estimates respectively the Strichartz 
estimates for wave equations with C 2 coefficients. The spatial localization scale 
is h'(\n(\x\))~2 \x\ while the frequency, instead of decaying, remains 0(1) due 
to the long range potential V. Hence the difference between P and its frozen 
coefficient version is 0(h'(\n(\x\))-$), which is more than the constant p 2 in 
the L 2 estimates. This is why we need also bounds on the second derivatives 
of a* 3 ', as required by pTheoremO 

Finally, the gradient potential can be added in as explained in the previous 
section. 

8. Appendix 
We consider a dyadic partition of unity in M. n , 

oo 

1 = S 

j=-oo 
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where Xj( x ) = Xoi^'^) is supported in \x\ ~ IK We also consider bump 
functions Xji x ) = Xji^^x) with slightly larger support, which equal 1 within 
the support of Xj such that XjXi = if \j — l\ > 2. 

Lemma 13. Let 1 < p < oo and —^<s<-. Then 



P 



\ u \\w 3 'P ~ ^ ] llXj^llii 



IP 



Proof. Let («.,•) be a sequence in W s,p . Arguing by duality it suffices to prove 
the bound 

oo 

ii xj u j\\w B -p ^ n^jiivi/ 3 ^ 

j=-oc 

With (D) = (1 + ID] 2 ) 1 / 2 we have 

oo 

II E XjUjWw^ = IKDy^XjUjllip 

j=-oo 

We write 



(d) s E = E • E c 1 -*)^)'* 

j=—oo j=—oo j=—oo 



Uj 



The terms in the first sum have almost disjoint supports and are easy to 
estimate. It remains to consider the second sum. We use bounds on the kernel 
of {D)~ s and its derivatives to estimate 

|(1 - Xj)(DYxMx)\ ^ IMIl,2 (s+ ^(2^ + \x\Y n - s . 
Then we conclude using 



E a,2^(2' + |*D— ||*,« E 



i ip n 
a,- r s > 

j=— OO j=— oo 



□ 



This is the main ingredient in the proof of 



Proposition 14. Let 5 > 0. Suppose that W G X (see Definition^). Then 
we have 

lim ||W||x(B(0,a)) = 

Proof. We assume that n > 3, the case n = 2 is similar. The result follows 
from the estimate 

(44) W\\x*\\XjW\\ ia p {X) 
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For one direction we write 

\(Wu,v)\ = \J2(XjWxjU,Xjv)\ 

% $^11x^11x11x^11 i ,2(^)11x^11 i 

* JJ r n + l ' n— 1 ^ n+1 ' n— 1 

~ II Xj^ll.g+l llXj^ll 2 (n + l) 1 2 (n + l) 2(n + l) 1 2(n+l) 

I 2 i n-l Wn+T' n-l i n-l W " + 1 ' 

< llXjV^H 2+1 ||U|| x 2(n+l) | |t>|| l 2(n+l) 

For the other, we consider separately sums with j even and with j odd: 



X! C*jW u j> u j) = Y (XjWxjUj,XjV, 

j even 

( W Y XjUj^Xj 



j even j even 



Vj) 



j even 



% w\w\\ Y xjUj\\,i,wi\\ Y WW i^^+ix 

' 4 \y n+l ' n—1 ' * p^n+1' n— 1 

jr e^en j even 

^5 X \\ U i\\ 2(n+l) 1 2{n+l) Kid 2(n + l) 1 2(n+l) 

/ n-l ^J/rT+T' n-l J n-l Wn+T' n-l 



□ 
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